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HYDRAULIC  MODEL  STUDY  OF  A  WATER 
INTAKE  UNDER  FRAZIL  ICE  CONDITIONS 

Thomas  ).  Tantillo 


INTRODUCTION 

The  Reynolds  Aluminum  Company  operates  a 
major  aluminum  reduction  plant  in  Massena,  New 
York.  The  facility  receives  its  cooling  water  through 
two  0.76-m-diam.  pipes  which  extend  107  m  into  the 
St.  Lawrence  River  T'n  se  pipe  inlets,  located  just 
downstream  of  the  Snell  Lock,  are  30  m  from  the 
southern  limit  ot  the  navigation  channel  (I  ig.  I).  Until 
a  stable  ice  cover  forms  on  the  river  upstream  of  the 
intakes,  they  are  subject  to  fra/il  icr  problems.  In 
the  past  tile  frazil  has  blocked  the  intake  entirely 
which  severely  curtailed  the  plant's  operation.  With 
the  onset  of  winter  navigation  the  stable  ice  cover  has 
not  been  allowed  to  form,  so  that  the  intake  has  been 
felt  t.)  be  subject  to  trequent  fra/il  problems. 

T ra/il  ice  occurs  when  turbulent  water,  initially 
above  0°C,  undergoes  cooling  at  a  rate  of  at  least 
0.01  °C/ 'hr  (Giffen  1973,  Granbois  1952).  As  the  rate 
of  water  cooling  increases,  larger  quantities  of  fra/il 
ice  arc  generated  because  of  the  greater  supercooling. 

The  most  common  method  of  keeping  trash  racks 
on  water  intakes  free  of  ice  accumulation  is  by  the 
application  of  heat.  There  is  much  disagreement  in 
the  literature  as  to  what  level  of  heat  is  necessary  per 
unit  area  of  trash  rack,  which  can  be  explained  by  the 
varying  water  velocities  and  bar  spacings  of  the  intakes. 
Rack  flushing  of  water  intakes  has  been  tried  at  a 
number  of  installations  with  little  success.  More  water 
is  used  to  back  flush  the  intakes  than  is  taken  in  be¬ 
fore  the  trash  racks  plug  again.  Other  installations 
use  traveling  screens  that  are  cleaned  automatically  or 
manually  (Giffen  19731. 

The  purpose  of  this  report  is  to  present  the  results 
ol  tests  on  two  protective  structures  that  were  modeled 
in  a  refrigerated  flume  to  minimize  the  fra/il  blockage 
at  the  intake.  The  refrigerated  flume  as  shown  in 
figure  2  is  described  in  Appendix  A. 


The  protective  structures  are  basically  concentric 
enclosures  around  the  intake,  which  arc  open  on  the 
top  and  which  set  up  a  vertical  velocity  profile  in  the 
water  entering  the  intake. 

MODELING  CRITERIA 

The  active  frazil  disks  produced  in  the  flume  are 
roughly  of  the  same  diameter  as  those  generated  in 
nature,  approximately  2-3  mm.  This  presents  a  pio- 
blcm  in  model  scaling  because  the  fra/il  disks  cannot 
be  scaled  down.  1  o  work  around  this  non-scaling  ot 
the  frazil  disks,  a  kinematic  model  is  proposed. 

The  kinematic  model  requires  that  the  model  and 
prototype  have  similar  velocities,  i.c. 

VJV  =  V  =  1  m 

p'  m  r  (II 

where  V  and  Vm  are  prototype  and  model  velocities 
respectively,  with  Vf  being  their  ratios.  A  suitable- 
model  scale  for  the  intake  structure  was  1 :24,  or 
length  ratio  Lr  =  24.  This  was  chosen  primarily  hv 
the  constraints  of  the  flume  width  and  water  depth. 
Since  the  internal  velocities  and  frazil  ice  particle- 
sizes  are  the  same  in  model  and  prototype,  the  drag 
forces  on  the  particles  after  they  enter  the  structure 
are  also  equal. 

Dimensional  analysis 

The  drag  force  on  a  frazil  particle  is  a  function 
of  the  water  velocity  V,  panicle  she  <i,  water  density 
p  and  water  viscosity  p: 

f  d  -  ( t',  4,p,  .p).  (2) 

Using  the  dimensional  equation 


Figure  I.  Loc  ation  of  Reynolds  Aluminum  Company  water  intake  on  St.  Lawrence 
River. 


Ld  =  K  VAdttpcp°  (3) 

eu  3  tan  be  dimensionallv  related  as 

F  =  (LF-])a  (L)B  (FT2L-4)c  (FTl.-2)0  (4) 

where  force  is  f,  length  L,  and  time  T. 

Solving  the  above  bv  use  of  the  n  theorem  results  in 

Fd  =  K  |  V2d2p(plpVd)°\  (5) 

This  yields  FJpV2d2  =  f(R()  where  the 
Reynolds  number  Re  is  equal  t opVd/p.  Thisf(Rf) 
can  be  defined  at  the  drag  coefficient  Crf,  and  White 
( 1 974)  gives  an  empirical  equation  for  the  drag  co¬ 
efficient  on  a  spherical  particle  as 

Cd  24/Rc  +  6/(1  +v/Rp  +  0.4.  (6) 


times  that  of  a  sphere  with  the  same  diameter  and 
low  Rc. 

A  computer  program  was  written  to  determine  the 
downward  vertical  velocity  at  which  frazil  disks  be¬ 
come  neutrally  buoyant.  This  velocity  will  occur  when 
the  drag  force  on  the  disk  equals  the  buoyant  force  on 
the  disk.  If  the  downward  velocity  component  is  low 
enough,  the  buoyant  force  will  overcome  the  drag 
force  and  the  disk  will  rise.  The  variables  are  the 
downward  vertical  water  velocity  frazil  disk  di¬ 
ameter  d,  and  thickness  of  frazil  disk  (typically  c//20). 
The  results  are  shown  in  Figure  3.  For  a  drag  force 
to  buoyant  force  ratio  (Fd/Fb)  less  than  1 ,  the  disk 
will  rise. 

For  frazil  disks  of  3  mm,  a  limiting  vertical  velocity 
of  1 .2  cm/s  is  reached  and  any  higher  velocity  will 
result  in  frazil  being  entrained  into  the  flow  and  into 
the  intake. 

I  rom  eq  1  the  water  discharge  scale  through  the 
intakes  Qf  is  thus 

<?r  =  z  V 


Work  done  by  (Hecn  (White  1974)  shows  that  a  disk 
normal  to  the  free  stream  has  a  drag  coefficient  0.85 


2 


(7) 


figure  3.  Drag  force /buoyant  force  ratio  vs.  disk 
diameter. 

with  the  time  ratio  T  merely 

T r  --  Lrl  Vr  =  Lr  (8) 

since  Vf  -  1  by  definition. 

F  low  patterns  around  the  structure  are  primarily 
the  result  of  inertial  and  gravitational  forces.  The 
F  roude  number,  which  indicates  the  relative  magni¬ 
tude  of  the  inertial  to  gravitational  forces,  is  there¬ 
fore  the  dimensionless  parameter  to  be  used  for  the 
free  stream  flow.  The  I  roude  number  is  defined  as 

Fm  -  Vls/gl.  (9) 

where  g  is  the  gravitational  constant.  The  F  roude 
numbers  of  the  model  Fm  and  prototype  Fp  should 
be  equal  to  ensure  similarity  for  the  free  stream  flow. 
The  free  stream  velocity  criterion  for  similarity  there¬ 
fore  is 


Vr  =  JTr 


(10) 


The  Saint  Lawrence  Seaway  Development  Corpo¬ 
ration  (Adams  1977)  obtained  current  direction  and 
velocities  by  using  drogues  in  the  South  Cornwall 
Channel  in  the  vicinity  of  the  Reynolds  Aluminum 
water  intake.  The  drogue  measurements  were  made 
on  1 5  and  20  )  uly  1977  with  flows  of  7050  and  7160 
m3/s  respectively.  The  results  show  that  the  mean 
velocity  over  the  intake  is  between  0. 1  2  and  0.30  m/s. 
The  average  value  of  0.21  m/s  was  used  for  this  study. 
The  model  free  stream  velocity  was  therefore  0.04  m/s 
for  a  1  :24  scale. 

The  model  and  prototype  scaling  parameters  are 
shown  in  Table  1  using  the  two  scaling  criteria. 


Table  1.  Scaling  results 


Prototype / 


Model 

Prototype 

Model 

Intake  flow 

245.8  cm 3 is 

0.142  m’/s 

-> 

('-r>‘ 

1 

Intake  velocity 
(cm/s) 

6.06 

6.06 

Frazil  disks 
(mm) 

2-3 

2-3 

1 

Intake  accumulation 
time 

1 

24 

'r 

Free  stream  velocity 
(m/s) 

0.04  3 

0.21 

Structural  scale 

1 

24 

/ 

Model  of  existing  intake 

A  1 : 24  geometric  scale  model  of  the  water  intake 
was  constructed  using  copper  tubing  for  the  inlet  pipe 
and  copper  sheet  for  the  conical  sections  and  baffles 
as  shown  in  Figure  4.  The  model  was  constructed 
from  prototype  drawings  supplied  by  Reynolds  Alu¬ 
minum  .  The  1 .83-m-diam.  prototype  intake  is  de¬ 
picted  in  I  igurc  5.  The  model  was  placed  in  the  cen¬ 
ter  of  the  flume  approximately  2  m  from  the  tailbox. 

The  model  intake  piping  system  was  calibrated  for 
flow  vs  differential  pressure  head  in  the  piping.  Man¬ 
ual  observation  of  a  manometer  recorded  the  changes 
in  flow  with  time. 


TEST  PROCEDURES 

The  scaled  intake  was  placed  in  the  flume  (see  Fig. 
6),  the  scaled  intake  (low  (?0  was  set  at  245.8  cm3/s. 
and  a  water  depth  of  35.6  cm  was  established.  Frazil 
ice  generation  was  initiated  and  time  vs  intake  water 
flow  was  recorded  manually  (Tantillo  1979).  This 
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Figure  4.  Geometric  scale  mode /  of  intake. 


Figure  5.  Section  ot  intake  structure. 


procedure  was  repeated  with  two  different  intake 
protective  structures  (one  circular  and  the  other 
square)  and  each  case  was  repealed  at  least  five  times. 

I  he  test  duration  was  limited  to  one  hour  due  to 
Ira/il  plugging  in  the  flume  tailbox. 

With  the  velocities  equal  in  the  model  and  proto¬ 
tv  pe  structures,  the  prototype  Reynolds  number  is 
24  times  as  large  as  that  of  the  model.  The  develop¬ 
ing  velocity  profile  would  therefore  be  different.  How¬ 
ever,  since  the  protective  structure  is  less  than  half 
a  pipe  diameter  in  length,  the  profile  in  both  cases  is 
lar  from  being  fully  developed  with  the  majority  of 
the  structure  area  having  the  core  velocity  with  a  small 
boundary  layer  near  the  diameter.  The  result  would 
be  a  slightly  lower  centerline  velocity  in  the  proto¬ 
tv  pe  than  in  the  model,  resulting  from  a  decreased 
boundary  layer  thickness  in  the  prototype  due  to  the 
higher  Reynolds  number  (Rohsenow  and  Choi  1961, 
Schlicting  1 955). 

Using  the  I  roude  criterion  for  scaling  the  free 
stream  velocity,  the  RL,  of  the  model  is  6200  whereas 
that  of  the  prototype  is  960,000.  f  or  the  round  pro¬ 
tective  structure  flow  separation  in  the  prototype 
would  take  place  farther  downstream  than  in  the 
model.  I  or  the  square  structure  the  separation  would 
occur  similarly  at  the  corners,  fhe  results  of  this 
would  be  less  of  a  wake  in  the  prototype  round  struc¬ 
ture  due  to  reduced  flow  separation,  and  similar  wake 
geometry  for  the  square  structure  in  the  model  and 
the  prototype  (White  1974,  Schlicting  1955). 


TEST  RESULTS 

The  unprotected  intake  structure  became  plugged 
quite  rapidly  with  fra/il  as  indicated  in  f  igure  7.  I  he 
average  time  was  1  7.2  min  or  6  hr  53  min  prototype 
time.  As  the  accumulation  of  ice  increases  in  the  in¬ 
take,  the  intake  flow  decreases  due  to  the  increasing 
head  loss  from  the  blockage. 

A  34.3-cm-diam.  x  20.3-cm-high  circular  protec¬ 
tive  structure  which  was  open  on  the  top  was  placed 
concentrically  around  the  intake  (fig.  8).  Using  the 
open  cross-sectional  area  the  initial  average  downward 
velocity  was  calculated  to  be  0.0027  m/s.  This  struc¬ 
ture  decreased  the  accumulation  rate  as  shown  in 
f  igure  9.  It  was  observed,  however,  that  this  config¬ 
uration  allowed  a  recirculating  vertical  eddy  to  de¬ 
velop.  This  nonuniform  flow  field  which  developed 
in  the  structure  allowed  fra/il  to  be  entrained  in  the 
higher  velocity  regions,  bringing  it  down  to  the  intake 
and  allowing  accumulation  to  occur. 

A  30.5-cni-square  structure  (20.3  cm  high)  having 
a  5-cm  45°  lip  was  tested  (I  ig.  10)  and  gave  an  initial 
average  downward  velocity  of  0.0027  m/s.  An  eddy 
still  developed  but  it  was  reduced  in  intensity  so  that 
lower  quantities  of  fra/il  were  entrained.  The  rate 
of  accumulation  on  the  intake  for  this  structure  is 
indicated  in  f  igure  1  1.  After  30  min,  which  corres¬ 
ponds  to  12  hours  for  the  prototype,  the  amount  of 
accumulation  is  reduced  by  an  average  of  50%  from 
that  of  the  round  protective  structuie.  This  is  due  to 
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f  igure  1 1 .  Flow  vs.  time  for  intake  with  square  protective  structure. 
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iQtirc  /  i.  -1<  cumulation  of  tra/il  it  p  after  25  minutes  on  unprotected  intake. 
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the  reduction  ol  the  eddy  in  the  protective  structure 
from  the  angled  lip. 

Typical  accumulations  of  tra/il  ice  are  shown  photo¬ 
graphically  in  I  igures  I  2  and  13.  Comparison  ot 
I  igures  I  2  and  I  3  shows  the  advantages  of  a  protective- 
structure. 


TEST  LIMITATIONS 

W  ith  the  present  Hume  configuration,  it  was  not 
possible  to  obtain  the  correct  I  roude  scaled  depth  and 
velocities  and  still  gener  ate  Ira/il  ice.  I  he  maximum 
depth  at  which  Ira/il  could  be  generated  at  the  correctly 
scared  stream  How  velocity  was  0. 875  nt  the  calculated 
sc. lied  depth.  I  hie  resulted  in  larger  recirculation 
patterns  in  the  model  protective  structures  than  would 
he  expected  in  the  prototype,  which  enhanced  the 
accumulation  ot  ha/il  ice  in  the  model,  (his  makes  the 
results  c  i  nisei  v  alive.  I  he  circulation  patter  ns  wci  e 
checked  at  the  o'rect  scaled  depth  without  Ira/it 
generation  anti  it  was  observed  that  there  were  lewet 
rec  itvi.'ation  eddies  at  this  greater  depth.  I  h is  would 
reduce  the  amount  ot  tra/il  being  drawn  to  the  intake. 

With  any  water  intake  protective  structure,  there 
must  he  a  minimum  amount  of  recirculation  in  the 
downward  direction.  Recirculation  will  pull  fra/il 
ice  into  the  intake  and  ultimately  will  result  in  block¬ 
age. 

The  modelingof  tra/il  ice  interactions  with  struc¬ 
tures  is  still  in  its  infancy.  There  is  much  yet  to  be 
leal  lied. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A  protective  structure  for  a  river  bottom  intake 


was  designed  which  reduced  the  vertical  water  velo¬ 
city  to  the  intake  enough  so  that  the  drag  force  of 
the  entrained  water  How  on  the  Ira/il  disk  was 
overcome  bv  the  buoyant  force  ol  the  tra/il  disk  • 

A  protective  structure  over  the  intake  which  cleated 
low  vertical  velocities  was  found  to  significantly  in¬ 
crease  the  time  before  the  intake  plugged  with  tra/il. 

It  is  recommended  that  studies  he  carried  one 
step  further  to  arrive  at  a  means  ol  reducing  Ira/il 
accumulation  h\  using  the  waste  heat  in  the  dis¬ 
charge  side  of  the  water  loop.  I  his  could  be  a 
more  cost-etlectivc  and  more  reliable  means  ot  cur¬ 
tailing  Ira/il  accumulation. 
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APPENDIX  A:  FLUME  TEST  FACILITY 

The  flume  facility  at  CRRFL  (Fig.  2)  is  36.6  m 
long,  1.2  m  wide  and  0.6  m  deep.  The  flume  has  a 
slope  range  of  -0.009  to  0.018  m/m.  The  system  has 
two  pumps  that  can  circulate  a  total  of  0.4  m  Is,  re¬ 
sulting  in  a  velocity  of  0.7  m/s  at  a  depth  of  0.46  m. 
This  flow  is  measured  with  inline  magnetic  flow 
meters  connected  to  flow  controllers.  The  valves 
allow  for  automatic  or  manual  flow  control. 

The  underside  of  the  flume  hottom  contains  two 


refrigerated  glycol  loops,  in  seven-  and  six- 
channel  paths,  through  which  opposing  or  parallel 
glycol  flows  can  he  run  in  either  direction.  Ihe 
glycol  temperature  can  be  adjusted  between  -18°C 
and  27°C. 

The  flume  is  situated  in  a  room  that  can  be  re¬ 
frigerated  to  -29° C.  It  has  two  /ones  of  cooling, 
each  with  four  air  handling  units,  which  are  fed  by  a 
liquid  ammonia  refrigeration  system.  A  refrigerated 
coil  in  the  water  storage  tank  can  chill  the  test  water 
to  0.3°C  prior  to  testing. 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Tantillo,  Thomas  J. 

Hydraulic  model  study  of  a  water  intake  under  frazil 
ice  conditions  /  by  Thomas  J.  Tantillo.  Hanover,  N.H. 
U.S.  Cold  Regions  Research  and  Engineering  Laboratory; 
Springfield,  Va.:  available  from  National  Technical 
Information  Service,  1981. 

ii,  12  p.,  illus.;  28  cm.  (  CRREL  Report  81-3.  ) 
Prepared  for  Directorate  of  Civil  Works,  Office  of 
the  Chief  of  Engineers/by  Corps  of  Engineers,  U.S. 

Army  Cold  Regions  Research  and  Engineering  Laboratory 
under  CWIS  31393. 
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